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I. INTRODUCTION 

The presence of hard multiple parton interactions in 
high-energy hadron-hadron collisions has been convinc- 
ingly demonstrated by the AFS [1], UA2 @, CDF 
and DO [5[ Collaborations, using events with the four- 
jets and 7 + 3-jets final states, thus, providing new and 
complementary information on the proton structure. 

In general, four high-i^T jets (or 7 + 3-jets) may be 
produced either in the collision of one pair of partons in 
2—^4 hard subprocess or via the simultaneous interac- 
tion of two parton pairs, that is in two 2—^2 subpro- 
cesses. In the last case in each dijet (or 7+ jet) system 
the transverse momenta of two high-Er jets more or less 
balance each other. 

The possibility of observing two separate hard colli- 
sions was proposed long ago. Their theoretical investiga- 
tion has a long history and goes back to the early days 
of the parton model [6-8] with subsequent extension to 
perturbative QCD [§^U 

Nevertheless, the phenomenology of multiple parton 
interactions relies on the models that are physically in- 
tuitive but involve significant simplifying assumptions. 
Therefore, it is extremely important to combine theo- 
retical efforts in order to achieve a better description 
of multiple interactions, in particular, double scattering, 
which will be the dominant multiple scattering mode at 
the LHC. In this letter we consider some steps towards 
this purpose. The cross section formulae currently used 
to calculate double scattering processes are revised in 
terms of the modified collinear two-parton distributions 
extracted from deep inelastic scattering (DIS). 

Let us recall, that with only the assumption of factor- 
ization of the two hard parton processes A and the 
inclusive cross section of a double parton scattering pro- 
cess in a hadron collision is written in the following form 



T (A,B) 



7T J r u(^i^2;bi,b 2 ;(3i,(3i) 



i,j,k,l 



x °tk ( x i > x 'i > Qi)°fi O2 ,x' 2 ,Ql) 
xT k i(x 1 ,x 2 ;b 1 - b,b 2 - h\Q\,Q\) 
xdxidx2dx 1 dx 2 d 2 bid 2 b2d 2 b, 



where b is the usual impact parameter, that 
is the distance between centres of incoming (the 
beam and the target) hadrons in transverse plane. 
Tij(xi, X2\ bi, b 2 ; Q 2 , Q 2 ) are the double parton distri- 
bution functions, depending on the longitudinal momen- 
tum fractions x\ and X2 and on the transverse position bi 
and b 2 of the two parton undergoing the hard processes 
A and B at the scales Qi and Q2. of k and a? t are the 
parton- level subprocess cross sections. The factor ra/2 
is a consequence of the symmetry of the expression for 
interchanging parton species i and j. m = 1 if A = B 
and m = 2 otherwise. 

The double parton distribution functions 
Tij(xi, X2\ bi, b 2 ; Q\ , Q|) are the main reason of 
interest in multiple parton interactions. These distribu- 
tions contain in fact all the information of probing the 
hadron in two different points simultaneously, through 
the hard processes A and B. 

It is typically taken that the double parton distribution 
functions may be decomposed in terms of longitudinal 
and transverse components as follows: 



r^(xi,x 2 ;bi,b 2 ;Qi,Q2) 
= D^(x 1 ,X2;QlQ 2 2)f(b 1 )f(b 2 ), 



(2) 



where /(bi) is supposed to be an universal function for 
all kind of partons with the fixed normalization, 



/(bi)/(bx - b)d 2 b l( l 2 b = / T(b)d 2 b = 1, (3) 



and T(b) = J /(bi)/(bi — h)d 2 bi is the overlap func- 
tion. 

If one makes the further assumptions that the longitu- 
dinal components D l £ (x±, x 2 ; Qi, Q2) re duce to the prod- 
uct of two independent one parton distributions, 

D^{x l ,x 2] Q\,Ql) = Dl{x 1 -Q 2 )D j h {x 2] Ql), (4) 

the cross section of double parton scattering can be ex- 
pressed in the simple form 



ra a (A) cr (B) 

2 CT ef f 



(1) 



°{A,B) 
CT eff = [Jd 2 b(T(h)) 2 



(5) 
(6) 
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In this representation and at the factorization of longi- 
tudinal and transverse components, the inclusive cross 
section of single hard scattering is written as 

i,k 

xTjfc^; t>i — b; Q^dxidx^bicPb 

i,k 

Y.D%{x 1 -Q 2 1 )f{b 1 - Wjdxxdx'^bidH 
= J2 J Diix^Q^a^ixux'^Dlix'^Q^dx^. (7) 

i,k 

These simplifying assumptions, even though rather 
common in the literature and quite convenient from a 
computational point of view, are not enough justified and 
should be revised [26|, |42|, |43| , while the starting cross sec- 
tion formula (pQ) was found (derived) in the momentum 
representation using light-cone variables in a number of 
works (see, e.g., Refs. 0, [H, HE US El) at the same 
approximations as in processes with a single hard scat- 
tering. 

II. MOMENTUM REPRESENTATION 

All the previous formulae were written in the mixed 
(momentum and coordinate) representation. Recall that 
in general for the case of the multiple parton interactions 
we have to use the Generalized Parton Distribution Func- 
tions (GPDF). In other words in the Feynman diagram 
(ladder) which describes the GPDF the parton momenta 
k>L (in the left part of diagram corresponding to the am- 
plitude A*) and Ur (in the right part of the diagram 
corresonding to amplitude A) may be different. Let us 
denote kh = k + q/2 and kn = k — q/2 where q is the 
momentum transfered through the whole ladder. Since 
the ladders in Figs. 1 and 2 form a loop we will call q — 
the loop momentum. In the previous formulae instead 
of transverse momentum q t we had used the conjugated 
coordinate b |62| . 

For our further goal the momentum representation is 
more convenient: 

a (A,B) = y / r^(xi,x 2 ;q;Qi,Q2) 

i,j,k,l J 

xaf k (x u x 1 )af l (x 2 ,x 2 ) 
d 2 q 

xT k i(x 1 ,x 2 ]-ci;Ql,Ql)dx 1 dx 2 dx 1 dx 2 -——r . (8) 

(2tt) 2 

The hard subprocesses A and B are originated by two 
different branches of the parton cascade. Note that only 
the sum of the parton momenta (in both branches) is 
conserved, while in each individual branch there may be 
some difference, of transverse (parton) momenta in 
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FIG. 1: A graph for double parton scattering due to the first 
term. A and B are the hard parton subprocesses. q is the 
momemtum transfered through the ladders L1,L2,L1 / and 
L2'. 
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FIG. 2: A graph for double parton scattering due to the 
second evolution term. A and B are the hard parton subpro- 
cesses. q is the momemtum transfered through the ladders 
L1,L2,L1 / and L2' 



the initial state wave function and the conjugated wave 
function. 

The main problem is in the correct calculation of 
these two-parton functions Tij(xi, x 2 \ q; Qf, Q 2 ) with- 
out the simplifying assumptions (|2]) and The 
case is that they are known in the current litera- 
ture [44]-[48| only at q = in the collinear ap- 
proximation. In this approximation the two-parton 
distribution functions, Tij(xi,x 2 ;q = 0;Q 2 ,Q 2 ) = 
Dfr(xi, x 2 ; Q 2 , Q 2 ) with the two hard scales set 
equal, satisfy the generalized Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations, derived 
for the first time in Refs. [13, EH, as well as single dis- 
tributions satisfy more known and cited DGLAP equa- 
tions [49l-[52|. The functions in question have a specific 
interpretation in the leading logarithm approximation of 
perturbative QCD: they are inclusive probabilities that 
in a hadron h one finds two bare partons of types i and j 
with the given longitudinal momentum fractions x\ and 
x 2 . 

The evolution equation for Tij contains two terms. The 
first term describes independent (simultaneous) evolution 
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of two branches of parton cascade: one branch contains 
the parton x\ and another branch — the parton x 2 . The 
second term accounts for the possibility to split the one 
parton evolution (one branch k) into the two different 
branches, i and j. It contains the usual splitting func- 
tion Pk->ij(z). The solutions of the generalized DGLAP 
evolution equations with the given initial conditions at 
the reference scales /i 2 may be written [24|, [29|, |53| in the 
form: 

= Dl\ h (x 1 ,x 2 ;^,QlQ 2 2 )+D^(x 1 ,x 2 ; f i 2 ,QlQ 2 2 ), 



Di\ h (x 1 ,x 2 ; fJ 2 ,QlQ 2 2 ) 



x^(^;m 2 ,Q?)^(?^ 2 ,Q1), 
D j ^{x u x 2 ^\QlQl) 



z\ 



Z 2 



min(Q?,Q2) 



1 — X2 



1-21 



= E 

hh'h' 



^ I d f s 



1 



xD J h (z! + z 2 ;v 2 ,k 2 ) — : — Pn> 



Z\ 



Z\ + z 2 
r x 1 



h^h'h' 



Z\ + z 2 



x 2 



xD£,(^-k\Qi)Dp(^k 2 : Qi) : 
J1 Z\ JZ z 2 

(9) 

where a s (k 2 ) is the QCD coupling, D 3 ^,{z \ k 2 , Q 2 ) are the 
known single distribution functions (the Green's func- 
tions) at the parton level with the specific £-like ini- 
tial conditions at Q 2 = k 2 . D^^ 2 (zi, z 2 , /j 2 ) is the ini- 
tial (input) two-parton distribution at a relatively low 
scale ii. The one parton distribution (before the slit- 
ting into the two branches at some scale k 2 ) is given by 

D(( Zl +z 2 ^ 2 ,k 2 ). 

Note that in Eq. (j9j) we assume that the loop mo- 
mentum q < fi is small and due to a strong ordering of 
parton transverse momenta during the collinear DGLAP 
evolution it may be neglected. 

The first term is the solution of homogeneous evolution 
equation (independent evolution of two branches) where 
in general the input two-parton distribution at low scale 
\i is not known. For this non-perturbative two-parton 
function at low z\ , z 2 one may assume the factorization 
D^ h \ Zl ,z 2 ,fx 2 ) ~ ( Zl ,fi 2 )D%\z 2 ,ii 2 ) neglecting 
the influence of momentum conservation {z\ + z 2 < 1). 

This leads to 



Dhifci, x 2 ]ii 2 ,Q\,Ql) 
Z^( W 2 ,Qf)I^( W 2 ,Q!). 



(10) 



Since the multiple interactions take place at relatively 
low transverse momenta, the contribution to the cross 



section from double scattering can be separated experi- 
mentally in the region of relatively low longitudinal mo- 
mentum fraction x\ and x 2 , where the factorization hy- 
pothesis (pTU]) for the first term is a good approximation. 
In this case the cross section for double parton scatter- 
ing can be estimated, using the two-gluon form factor of 
the nucleon F 2g (q) [26|, [54| for the dominant gluon-gluon 
scattering mode (or something similar for other parton 
scattering modes), 



D,lxl 



i,j,k,l ' 



. j Di{xv,H 2 ,Q\)D j h {x 2 ;ii 2 ,Ql) 
xaf k (x 1 ,x' 1 )af l (x 2 ,x' 2 )D^(x' 1 ;iJ, 2 ,Q 2 )D l h ,(x , 2 ;ii 2 ,Q 2 2 ) 



d 2 c 



(11) 

From the dipole fit F 2g (q) = l/(q 2 /m 2 + l) 2 to the two- 
gluon form factor follows that the characteristic value of 
q is of the order of gluon mass m g , and therefore the 
initial conditions for the single distributions can be fixed 
at some reference scale \i ~ m g because of the weak log- 
arithmic dependence of these distributions on the scale 
value. In this approach J F^g (q) gi yes the estima- 
tion of [cr e ff] _1 . 

The second term in Eq. (|9)) is the solution of complete 
evolution equation with the evolution from one "nonper- 
turbative" parton at the reference scale. Here the in- 
dependent evolution of two branches starts at scale k 2 
from a point-like parton j' . Therefore we have no the 
form factor F 2g (q) which suppresses the large q t domain. 
The corresponding contribution to the cross section reads 



^ J dx\dx 2 dx 1 dx 2 J 



D,2x2 
a (A,B) 

min(Q?,Q2) rf2 



i,j,k,l 



l-2i 



x E 

hh'h' 



2 a s (k 2 ) f dz x f dz 2 
2tt/c 2 J zi J z 2 



xD J h ( Zl +z 2 ;[i 2 ,k 2 ) 



z\ + z 2 



-Pi' 



Z\ 



j'^h'h' 



Z\ + Z 2 



x E 

3'h'h' 



xI^-;^Q?)^(-;* 2 ,Ql) 

J Z\ J2 z 2 

xaf k (x 1 ,x / 1 )af l (x 2 ,x 2 ) 
i-4 i-,i 



dk 



xDl^+z^v^k' 2 ) 



> 2 a s (k 2 ) 
2itk' 2 



f dzi f dz 2 
J Zi J z 2 



z\ + z 2 



Pj'^h'h' 



Z\ + z 2 



— \k , (Dn 



). (12) 



or in less transparent, but substantially shorter form: 

D,2x2 



a 



— J dxidx2dx 1 dx 2 

i,j,k,l 



(A,B) 



(2tt)2 



xD l ^ 2 (x 1 ,x 2 ]q 2 ,Q 2 1 ,Ql)af k (x 1 ,x 1 ) 



xaf l (x 2l x 2 )D k h t /2 (x , 1 ,x , 2 ;q 2 1 Q{, Q l 2 ). (13) 

There is yet the combined ("interference") contribu- 
tion, which is written by analogy, 



— / dx\dx 2 dx x dx 2 \ 



Fl(q) 



D,lx2 
T (A,B) 

d 2 q 
(2ir) 2 



x [D % h (xi; n 2 , Q\)D :i h (x 2 ; /j 2 , Ql)af k (x!, x' r ) 
x af t (x 2 ,x' 2 )D^ 2 (x' 1 ,x' 2 ;q 2 ,Q 2 1 ,Ql) 
+D% 2 (x 1 ,x 2 ;q 2 , Ql, Q\)at k {x\, Xi) 
xaf t {x 2 , x 2 )D\, (x[ ; /x 2 , Q\)D l h , (x' 2 ; /x 2 , Q\)]. 

(14) 

The equations ([TT]) . (p~3|) and (fT4|) are our solution of 
the problem — the estimation of the inclusive cross sec- 
tion for double parton scattering, taking into account the 
QCD evolution, in terms of the well-known collinear dis- 
tributions extracted from deep inelastic scattering. How- 
ever, here one should note that the input two-parton 

distribution D 3 ^' 32 (z±, z 2l /i 2 ) may be more complicated 
than that given by factorization ansatz ([TO]) . Let us dis- 
cuss in more detail the second term, that is the 2x2 
contribution. 



III. 



DISCUSSION AND CONCLUSIONS 



The contribution to the cross section from the second 
term induced by the QCD evolution does not reduce to 
the simple contribution with a some new, other, con- 
stant effective cross section <j e ff, as it was done in first 
estimations pi! E3, EH • The QCD evolution effects for 
the cross section are anticipated to be larger than for the 
two-parton distribution functions, for which they were es- 
timated in Refs. [11,113 on the level of 10 % -30% in com- 
parison with the "factorization" components at x ~ 0.1 
and Q ~ 100 GeV. Indeed, in Eq. (fl2]) the integration 
over q has no a strong suppression factor F 2g (q) and the 
phase space integral may be estimated as, 



Q 2 r Q' aiA 
dq 2 



d^_ r c 

k 2 L 



dk 



'2 



k' 2 



2Q 2 , Q 2 >/i 2 , (15) 



where within the leading order (LO) accuracy we take q 2 
as the lower limit in k 2 and k 2 integrations; at q 2 > k 2 



the loop momentum q t destroys the logarithmic structure 
of the integrals in collinear evolution from k 2 to Q 2 . 

We see that for a large final scale Q 2 the second (2 x 
2) contribution dominates being proportional to Q 2 in 
comparison with the 1 x 1 or 1 x 2 components which 
contributions ~ m 2 ~ l/a e ff are limited by the nucleon 
(hadron) form factor F 2g |63j . 

The real gain is, of course, less, since the running cou- 
pling constant and the distribution functions are loga- 
rithmically dependent on the integration variables and we 
have the additional suppression factor, which is inversely 
proportional [HH to the initial gluon and quark multiplic- 
ities squared; due to the different "normalization": the 
second term evolves from one "nonperturbative" parton 
unlike the first factorization term having the two initial 
independent "nonperturbative" partons at the reference 
scale. 

As a result, the experimental effective cross section, 
a eff P ' wn i cn is n °t measured directly but is extracted, 
using the normalization to the product of two single cross 
sections: 



a 



7 +3j 
DPS 



exp-i-l 
eff J 



(16) 



in both the CDF and DO experiments, should be depen- 
dent on the probing hard scale. It should decrease with 
the growth of the resolution scale due to the fact that 
all additional contributions to the cross section of dou- 
ble parton scattering are positive and increase. Here a 13 
and a JJ are the inclusive 7+ jet and dijets cross sections, 
a Dp 3 s ^ s the inclusive cross section of the 7 + 3 jets events 
produced in the double parton process. It is worth notic- 
ing that the CDF and DO Collaborations extract <r^ p 
without any theoretical predictions on the 7+ jet and di- 
jets cross sections, by comparing the number of observed 
double parton 7 + 3 jets events in one pp collision to 
the number of 7 + 3 jets events with hard interactions 
occurring in the case of two separate pp collisions. 

The recent DO measurement [5[ of this effective cross 
section, <r^ p , done as a function of the second (ordered 

in the transverse momentum pr) jet pt, Vt^ ' •> that 
can serve as a resolution scale, shows a tendency to 
be dependent on this scale. In Ref. [56] this fact was 
interpreted as the first indication to the QCD evolution 
of double parton distributions. 

We have to emphasize that the dominant contribu- 
tion to the phase space integral (fT5j) comes from a 
large q 2 ~ Q 2 and strictly speaking we have no place 
for the collinear (DGLAP) evolution of two indepen- 
dent branches of the parton cascade (i.e., in the lad- 
ders L1,L2,L1 / and L2') in the 2x2 term. Formally 
in the framework of collinear approach this contribution 
should be considered as that caused by the interaction 
of one pair of partons with the 2—^4 hard subpro- 
cess [64]. Recall however that in the estimate (fT5|) we 
neglect the anomalous dimension, 7, of the parton dis- 
tributions Dj(x/ 'z,k 2 \Q 2 ) ex (Q 2 /k 2 )^. In collinear ap- 
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poach the anomalous dimensions 7 ex a s « 1 are as- 
sumed to be small. On the other hand in a low x re- 
gion the value of anomalous dimension is enhanced by 
the \n(l/x) logarithm and may be rather large numeri- 
cally. So the integral over q 2 is slowely convergent and 
the major contribution to the cross section is expected 
to come actually from some characteristic intermediate 
region, m 2 « q 2 « Q 2 (Qi < Q2), with the not such 
strong sensitivity to the upper limit of ^-integration as in 
the case of the phase space integral. Therefore it makes 
sense to consider the numerical contribution of the 2x2 
term even within the collinear approach at the LHC kine- 
matics with the available large enough Q\ and Q2 with 
respect to m g (providing the wide enough integration re- 
gion and the characteristic intermediate momenta q). 

Next in a configuration with two quite different scales 
(say, Q\ « Q 2 ) the upper limit of q 2 integral is given 
by a smaller scale (at q > Qi the hard matrix element 
corresponding to a A starts to fall down with q t ). In this 
case the collinear evolution from the scale q = Qi up to 
the scale Q2 in the ladders (parton branches) L2 and 
L2' is well justified. 

Moreover the integrals over Z2 in Eq. ([T2]) are 
concentrated at a rather large z giving enough space 
for the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolu- 
tion [57M6QI] . Thus it may be intersting (and even more 



justified theoretically) to consider the multiple parton 
interactions in the framework of the BFKL approach 
(see the recent paper [6lj). 

In summary, we suggest a practical constructive me- 
thod how to estimate the inclusive cross section for dou- 
ble parton scattering, taking into account the QCD evo- 
lution, in terms of the well-known collinear distributions 
extracted from deep inelastic scattering. We support also 
the conclusion of Refs. [56|, [6l| that the experimentally 
measured effective cross section, cr^ p , at a normaliza- 
tion ([T6]) and the presence of the evolution (correlation) 
term in the two-parton distributions, should decrease 
with the growth of the resolution scale Q 2 . 
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